Transfer of BECs through discrete breathers in an optical lattice 
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We study the stability of a stationary discrete breather (DB) on a nonlinear trimer in the frame- 
work of the discrete nonlinear Schrodinger equation (DNLS). In previous theoretical investigations 
of the dynamics of Bose-Einstein condensates in leaking optical lattices, collisions between a DB 
and a lattice excitation, e.g. a moving breather (MB) or phonon, were studied. These collisions 
lead to the transmission of a fraction of the incident (atomic) norm of the MB through the DB, 
while the DB can be shifted in the direction of the incident lattice excitation. Here we show that 
there exists a total energy threshold of the trimer, above which the lattice excitation can trigger 
the destabilization of the DB and that this is the mechanism leading to the movement of the DB. 
Furthermore, we give an analytic estimate of upper bound to the norm that is transmitted through 
the DB. Our analysis explains the results of the earlier numerical studies and may help to clarify 
functional operations with BECs in optical lattices such as blocking and filtering coherent (atomic) 
beams. 

PACS numbers: 67.85.De, 63.20.Pw, 03.75.Lm, 42.65.Tg 



1 . Introduction 

Since the experimental realization of Bose-Einstein Con- 
densation (BEC) of ultra-cold atoms in optical lattices 
(OLs) [lj, many groups of researchers have achieved an 
extraordinary level of control over BECs in optical traps 
[2-5J. Among other important applications, this control 
has allowed for the investigation of analogues of com- 
plex solid state phenomena [6HTD]. Technologically, the 
emerging field of "atomtronics" promises a new genera- 
tion of nanoscale devices. 

An important generic feature of nonlinear lattices is the 
existence of discrete breathers (DBs), which are spatially 
localized, time-periodic, stable (or at least long-lived) ex- 
citations in spatially extended perfectly periodic discrete 
systems [TTHT3] . DBs arise intrinsically from the combi- 
nation of nonlinearity and the discreteness of the system. 
DBs have been observed in a variety of systems, includ- 
ing Josephson-j unction arrays [14, 15J, micromechanical 
systems [16] , nonlinear waveguide arrays [171 [18] , a he- 
lix proteins [19J, spins in antiferromagnetic solids [20|l2T] 
and BEC |6J. The existence, stability, and other prop- 
erties of DBs have been studied theoretically throughout 
the last decade [T3l [22H24] . Among other results, it was 
shown that they act as virtual bottlenecks which slow 
down the relaxation processes in generic nonlinear lat- 
tices (25H2BJ. 

Many theoretical studies of the dynamics of a BEC 
trapped in an OL use the discrete nonlinear Schrodinger 
equation (DNLS) to model the system. Several recent 
studies based on the DNLS have observed the collision 
of a stationary DB with a lattice excitation, e.g. a mov- 
ing breather or phonon [27-29J. If the amplitude of the 
lattice excitation is too small, it is reflected entirely from 
the DB. Above a specific threshold amplitude, part of 



the incident atomic norm is transmitted through the DB 
while the DB is destabilized and shifted by one or few 
lattice sites towards the incoming moving breather (MB) 
|48| . This transmission process plays a central role in 
the occurrence of scale-free atomic avalanches observed 
for a whole range of nonlinearities in leaking optical lat- 
tices [28J. However, this process has heretofore not been 
understood analytically. 

In this article we analyze, analytically and numerically, 
the collision process of a stationary DB with a lattice ex- 
citation. To this end, we study the nonlinear trimer, 
i.e. the DNLS with M = 3 lattice sites (see e.g. [5] for an 
experimental realization of a similar system). We calcu- 
late analytically the threshold for the destabilization of 
the DB as well as an upper bound to the atomic norm 
that can be transmitted through the DB. The threshold 
and the transmission process are described by introduc- 
ing a 'Peierls-Nabarro energy landscape' which restricts 
the accessible region of the phase space for excitations on 
the trimer. The 'local Ansatz', [30| l3T]. an approach ap- 
plied successfully to studies of DBs on nonlinear lattices, 
suggests that the results we find for the trimer will apply 
to extended lattices; the agreement of our results with 
the recent numerical studies [28J confirms this. For the 
rather large nonlinearity we shall consider in the sequel, 
DBs are well localized, and the most basic and important 
DBs occupy only three sites. Within the local Ansatz we 
consider only the central site and the two neighboring 
sites of a DB. This allows us to reduce the high dimen- 
sional dynamical problem involving M sites to the non- 
linear trimer. A detailed analysis of the reduced problem 
[3Q| l3T] shows that (i) the DB corresponds to a trajectory 
in the phase space of the full system which is practically 
embedded on a two-dimensional toroidal manifold, thus 
being quasiperiodic in time; (ii) the full DB can be re- 
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produced accurately within the nonlinear trimer approx- 
imation. 

Although we focus here on BECs, our results are also 
relevant in a wide range of other contexts in which the 
DNLS applies, most prominently coupled nonlinear opti- 
cal waveguides [171 \T$\ . 



2. The Model Hamiltonian 

The Bose-Hubbard Hamiltonian is arguably the sim- 
plest model that captures the dynamics of a dilute gas 
of bosonic atoms in a deep optical lattice, with chemical 
potential small compared to the vibrational level spac- 
ing (see e.g. |3J for a review). In the case of weak in- 
teratomic interactions (superfluid limit) and/or a large 
number of atoms per well (so that the total number of 
atoms N ~ (9(10 4 — 10 5 ) is much larger than the num- 
ber of wells M), a further simplification is available since 
the BECs dynamics admits a semi-classical (mean-field) 
description. The resulting Hamiltonian describing the 
dynamics is 



M 



M-l 



n = J2Wn\ 4 ^^n\ 2 ]- 1 Z ^(CVWl+C.C.), (1) 



where n = 1,...,M is the index of the lattice site, 
\il) n (t)\ 2 = N n (t) is the mean number of bosons at site n 
(also referred to as the norm N n (t)), U = Ai:h 2 a s V e ^/m 
describes the interaction between two atoms at a single 
site (14ff is the effective mode volume of each site, m is 
the atomic mass, and a s is the s-wave atomic scattering 
length), ji n is the on-site chemical potential, and J is 
the tunneling amplitude. The "wavefunctions" ifj n can be 
used as conjugate variables with respect to the Hamilto- 
nian %, leading to a set of canonical equations 



dr 



dH 



dH 



(2) 



which upon evaluation yields the Discrete Nonlinear 
Schrodinger Equation 



A|VVi| 2 VVi - \ [^n-l +^n+l] • (3) 



Here, A = 2/7/ J is the nonlinearity and t = Jr is the 
normalized time. In Eq. [3] we have set \i n = Vn. 

The DNLS can be applied to a remarkably large va- 
riety of systems, in particular this mathematical model 
describes (in the mean-field limit) the dynamics of a BEC 
in an OL of size M [36]. Experimentally, the tunneling 
rate J can be adjusted by the intensity of the standing 
laser wave field. A powerful tool to modify the on-site in- 
teraction U is via a Feshbach resonance, where the atomic 
interactions can be controlled over a large range simply 




Figure 1: DB solutions for the symmetric case ipi = ip3 
(Eq. |7])). For a nonlinearity A > 5.04 four symmetric so- 
lutions exist. The solution for iV~2 > 1/2 is termed a bright 
breather, while the solution for N2 — > corresponds to a dark 
breather (see Eq. (12)). The dashed vertical line at N2 = 1/3 
marks the asymptote for the phase-wise and antiphase-wise 
time-periodic solutions for A — » 00. 



by changing a magnetic field. A Feshbach resonance in- 
volves the coupling of free unbound atoms to a molecular 
state in which the atoms are tightly bound. When the 
energy levels of the molecular state and the state of free 
atoms come closer, the interaction between the free atoms 
increases. Thus, the nonlinearity A can be varied experi- 
mentally. Here we will treat the repulsive case explicitly 
(A > 0); however, the attractive case can be obtained via 
the 'staggering' transformation ip n — ^ ( — l n )tfj n followed 
by time-reversal t —t [13J. 



3. Equations for the Nonlinear Trimer and Asymptotic 
Solutions 

To analyze the transfer of norm through a DB (and 
related the stability of the DB) during a collision in the 
nonlinear trimer (M = 3) we begin with the equations 



id t i/j 2 = A|^2| 2 ^2 - ttOAi +^3) 



^2 



idt^z = A|^ 3 | ^3 - ^2 • 



(4) 



We normalize the wave functions such that the total 
atomic population reads 



A I 



N = J2\lPn\ 2 = l- 



71=1 



To find single frequency breather solutions in Eq. Q for 
the symmetric case ipi = ^3, we assume 



ip n (t) = A n e 



iwt 



(5) 
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with amplitudes A n and frequency w. This Ansatz, to- 
gether with the conservation of particle number, leads to 
the set of equations 



wAi = \A\ - -A 2 
-wA 2 XAl - Ax 



1 = 2A\ 



Ah 



(6) 



Let us first calculate the relation between the (atomic) 
norm N 2 = A 2 at the central site and the nonlinearity A. 
From Eq. do) , we find 



\(N 2 ) = ± 



V2(2N 2 



1) 



^N 2 (l-N 2 )(W 2 - 1) 



(7) 



We have four solutions above the bifurcation point at 
A « 5.04 and two solutions for < A < 5.04 (see Fig. [I}. 
To gain further insight into the nature of the symmetric 
solutions in the t rimer, we will revisit Eq. (|6|, which we 
convert into a quartic equation 



where 



-^=x 3 - V2\x 
V2 



cos(arctan(x)) 



1 = 0, 



sign(x) 



vT 



Expansion of the exact real solutions of Eq. 
the limiting case A — >• gives 

A 5 a 2 + 0(A 3 ) 



(8) 

(9) 
in A for 



X\ — 1 



x 2 = — 1 



4y/2 
A 



4V2 64 
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+ ^A 2 + 0(A 3 ). 



(10) 



At A = 0the solution ^(t) = (^i(t), ^ 2 W, ^3 (*)) of Eq. (|6j) 

at time t = reads ^ xliX2 )(0) = (1/2, ±1/^2, 1/2). 

The antisymmetric breather configuration ^(0) = 

(-1/^2,0,1/^2) is not included in our Ansatz, as we 

restrict ourselves to symmetric solutions. 

Expansion around the limit A — » oo leads to four real 

solutions 

1 1 1 \+o(\- r >) 



x 3 



X4 



X 5 



x 6 



V2\ 
1 
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1 
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V2 A 
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14v / 2^3+0(A- 5 ) 
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WV2 A 2 
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(11) 



2\/2A WV2X 2 
For infinite A the solutions of Eq. (|6]) at time t = are 

$*3)(0) = (0,1,0) 
^ (x4) (0) = ( v ^72,0,v / l72) 
^(x5,x6)(0) = (l/>/3,Tl/V3,l/V3), (12) 
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Figure 2: (a) The lower part of the PN energy landscape 
exhibits three minima separated by saddle points, (b) The 
phase space of the trimer is restricted to lie between the two 
parts of the PN shell, which consists of the lower and upper 
part of the PN landscape, which are shown in this panel, (c) 
Contour plot of (a), the lower part of the PN energy landscape. 
The three minima and saddle points are clearly visible. The 
minimum at A\ — A3 = 0.17 in (c) corresponds to the bright 
breather. The figure is plotted for A = 3. 



where the solution ^03) is called a bright breather, ^(^4) 
is named a dark breather (due to lack of norm at the 
central site) and ^05, ^6) are phase- wise and antiphase- 
wise time-periodic solutions. 



4- The Peierls-Nabarro energy landscape 

Having found the symmetric DB solutions, we next fo- 
cus on the transfer of norm through a bright breather, 
where the stability of the breather will play a crucial 
role. We start by introducing the concept of a 'Peierls- 
Nabarro (PN) energy landscape'. It is related to the PN 
potential, which reflects the fact that discreteness breaks 
the continuous translational invariance of a continuum 
model [371 EBJ • The amplitude of the PN potential may 
be seen as the minimum barrier which must be overcome 
to translate an object by one lattice site. Regarding DBs, 
the Peierls-Nabarro barrier is given by the energy differ- 
ence \E C — Eb\, where E c and E^ are the energies of a DB 
centered at a lattice site and between two lattice sites. 

We define the Peierls-Nabarro energy landscape as fol- 
lows: for a given configuration of amplitudes, A n , the 
PN energy landscape is obtained by extremizing H with 
respect to the phase differences 5<pij = — : 



Hp N = mm(-H) ; max (-if) , 



(13) 
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where ijj n = A n exp(z0 n ) and H pN and i^p N are the 
lower and upper part of the PN landscape. As we will 
see later, the bright breather solution ^03) 1S located at 
an extremum of Hp N . The minus sign in the definition 



(13) was added for convenience to assure that the bright 



breather is found in a minimum (and not in a maximum) 
of the lower PN landscape. The phase differences ex- 
tremizing the Hamiltonian 

H = ±(A$ + Al + Ai) 

-(A!A 2 cosOi -<j> 2 ) + A 2 A 3 cos(0 2 -0s))(14) 

are 5<pi 2 = 5<j>23 € {0, n}. Hence, the upper and the lower 
PN energy landscapes read 



A 



PN 



(A 4 1 + A A 2 + Ai) + (A 1+ A 3 )A 2 (15) 



and 



PN 



■ A 4 2 + At) - (A, + A 3 )A 2 . 



(16) 



In Fig. [2] the PN landscape is visualized for A = 3. 
The PN 'shell', consisting of the upper and lower land- 
scapes, bounds the phase space of the trimer. Since the 
DB whose properties we are studying corresponds to a 
minimum on i^™, we shall focus on this landscape. As 
shown in Fig. [2] (c), the projection onto the A\ — As 
plane exhibits three minima which are separated by sad- 
dle points (called 'rims' in the following). For A — >• oo 
the saddle points are located at A\ = A2 = <\/l/2 (which 
in the following will be the saddle point of interest) and 
A2 = As = a/1/2. The energy threshold E t h rs at the rim 
(obtained from Eq. (16)) reads 



The following investigation holds for an effective non- 
linearity in a range around A = X/M ~ 1 [49J, which is in 
the critical regime where scale-free avalanches of BECs 
were found in [28j. 



5. The Threshold for transfer of norm 

To study the influence of the PN landscape on the 
stability and the transfer of atoms through the DB, we 
first consider the fixed point corresponding to the bright 
breather. The initial amplitudes A\ are obtained by in- 
verting Eq. ([7]) for N2 > 1/2; hence an initial condition 
for the bright breather reads ^ 6 (0) = (-A\, A b 2 , -A\). 
Then perturbations are added to site 1. In dynamical 
systems terminology, the phase space of the trimer is 
'mixed', consisting of regular islands surrounded by the 
chaotic sea [28, 39J. DBs are located inside the regular 



islands of the phase space, provided that their frequency 
(and multiples of their frequency) lie outside the phonon 
spectrum [HU| l3Tj . If a perturbation is large enough, 
it can push the orbit out of the regular island into the 
chaotic sea, destabilizing the DB. 

We now use the following initial condition: 



W) = {-{A\ + 8 A )e iS *,A 2 ,-A\), 



(18) 



where A2 = (1 — — 1^3 1 ) ' ensures total norm 

N = 1. Compared to the bright breather, we have added 
an amplitude 5 a to site 1 and the phase is rotated by 
8(f). The initial condition (18) is visualized in Fig. fright 
panels), where we have fixed 5^ = tt and increased 5 a 
in (a-d). Note that although the phase rotation does 
not alter the norms |^i| 2 ? it drastically changes the total 
energy of the trimer which we define as E t = —H (see 



Eq. g4j). 

In Fig. j3^1eft) we show the dynamics for increasing 
total energy E t , where the arrows on the orbits (black 
curves) mark the direction of time. For E t < £?thrs the 
areas in phase space are disconnected, leading to sub- 
threshold dynamics depicted in Fig. [3^. 

In contrast, for E t > £ahrs the orbit is allowed to pass 
the rim such that the majority of the norm migrates from 
site 2 to site 1 (Fig.^), while norm is transferred to site 
3. The larger E t (Fig. [3]:), the larger is the size of the 
"bubble" (by the term bubble we denote the accessible 
region of the PN landscape for A\ > l/y/2). Hence, as 
we see from Fig. [3]>c(left), an upper limit to the norm 
that can possibly be transmitted through the DB can be 
read from the maximum value of As inside the bubble 



4* 
^3 



max As 

At>l/V2 



(19) 



for fixed total energy E t . 

Finally, for even larger E tl the orbit visits large parts 
of the phase space and large amplitudes Ai(t) are found 
at all three sites, as depicted in Fig. where most of 
the orbit resides in the chaotic regime (Fig. |3]i). In this 
case, no controlled shift of the DB from site 2 to 1 is ob- 
served and the DB becomes dynamically unstable. These 
dynamical instabilities can be associated with a (partial) 
depletion of the BEC [40, 41J; a detailed study of these 
effects is beyond the scope of our present investigation. 

The upper bound A3 can be calculated analytically 
noting that for the projection onto the A\— As plane, con- 
dition dAs/dAi 
leads to 



holds. Implicit derivation of Eq. 16 



= 1- 2A\ - A x As 



■Ai 



+2AAia/1-A?-A§(-1 + 2A? + j4§), (20) 



that determines the maximum value of As in the "bubble" 
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Figure 3: Dynamics on the PN landscape for increasing total energy of the trimer. (a) (left) A contour plot of the lower PN 
energy landscape Hp N is shown for total energy below the rim (E t — —1.32 < i^thrs = —1.311). A projection of the orbit 
onto the A1-A3 plane is over-plotted (black curve), (middle) The corresponding amplitudes Ai(t) indicate that the maximum 
amplitude remains at the central site. The dashed vertical line marks the time interval [0, 25] for which the orbits in the left 
picture are plotted, (right) A sketch of the initial condition shows that the excitation at site 1 is slightly below threshold, 
(b) Destabilization of the DB for Et = —1.310 > E t hrs just above the rim. We see that the rim of the PN landscape clearly 
restricts the dynamics and governs the destabilization process of the DB. (c) For higher total energy E t — —1.28 the bottleneck 
at the rim widens and the maximum norm transmitted to site 3 is increased, (d) For even higher total energy Et — —1.04 the 
orbit explores large parts of the phase space and visits all three sites. The grey shaded areas are forbidden by the upper PN 
landscape Hp N . In all cases A = 3, 5$ = 7r. For other values of 5$ the same qualitative behavior is found. 
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Figure 4: Norm transfer through a bright breather. The max- 
imum (atomic) norm at site 3 detected after the collision of 
the DB with a lattice excitation is shown as a function of 
the total energy of the trimer. The discrete symbols indicate 
three different values of the initial phases. For Et < E t hrs the 
DB is stable and practically no transfer of norm takes place 
on short timescales. For E t > ^thrs we observe instability of 
the DB centered at site 2: The breather migrates to site 1 and 
norm is transferred to site 3. An upper bound to max(A^3(t)) 
is calculated from the PN landscape, both analytically (dot- 
ted dashed line, cf. Eq. (21)) and numerically (solid line). 
The analytical calculation is performed in the limit for large 
A and therefore deviates slightly from the exact numerical re- 
sult. We used 5023 = tt which is the typical case observed in 
[28| for DBs in an extended leaking optical lattice and A = 3. 



of the PN landscape as 



A* 3 (S,X) = 



1 

i 



V(l - 25) (1 - 28) 
2V\ 4V2X 
2(5 -5 2 ) 25 - 20(<5 



5 2 ) 



4VT- 

11/8- 



25X 3 / 2 
-35-S 2 



64A 2 
7(<5 3 - 7<5 4 ) 



2(1 - 2(5) 3 / 2 A 5 / 2 



+o(\- 



(21) 



where 5 = E t — E t h rs > is the energy relative to the 
destabilization threshold. Without the A -5 / 2 term the 
exact value is underestimated. If we truncate the ex- 
pression after the A -3 / 2 term, the deviation from the 
exact result roughly doubles compared to what is shown 
in Fig. [4j 

How general is the transfer mechanism that we de- 
scribe? In Fig. [4j the maximum norm N3 that is found 
at site 3 after the transfer of atoms through the DB is 
shown as a function of the total energy for three ini- 
tial phase differences = tt — 5$ = 0, 7r/4, tt/2. The 
dashed vertical line at E t = E t hrs — —1.311 marks the 
total energy at the rim and is identified with the desta- 
bilization threshold of the DB. Evidently, the transfer 
mechanism found does not depend on parameters 5 a and 
5(f) individually, but rather on the total energy which 
determines the accessible region of the PN landscape. 
Moreover, the transfer mechanism itself appears nearly 



independent of the choice of the initial phase difference 
5(j)\2. The maximum norm detected at site 3 is closely be- 



21 
1.2 



low the upper bound Ng(5) = 1^45 (<5) | 2 given by Eq, 
which holds for E t < —1.2. For increasing E t > - 
the fluctuations of max(A^) for different orbits with sim- 
ilar total energy become larger, because the orbits ex- 
plore a larger part of the phase space. As a consequence, 
max£<T(A/3(£)) depends on the chosen time interval [0, T] 
in which the maximum norm is detected (In all cases we 
set T = 100). In contrast, for lower total energy where 
a nontrivial upper bound N£(5) of the transferred norm 
holds, max^<T(A^3(t)) does barely depend on T already 
after very few oscillations of A 3 (t) (cf. Fig. |3b-c( middle)). 



6. Connection to PN barrier 

In order to gain further insight into the relation be- 
tween the rim of the PN energy landscape and the PN 
barrier, let us consider a nonlinear trimer where we omit 
the nonlinear on-site interaction term at site 3. The equa- 
tions of motion read 



A|^ 2 | 2 ^2 - -(V^i +^3) 



(22) 



Using the initial condition (18) we find that on the 
timescale where the destabilization process of the DB 
centered at site 2 takes place, the dynamics is not 
changed substantially compared to the results for the 
nonlinear trimer (see Fig. [5|. Hence, in order to describe 
the destabilization of the DB (and the basic mechanism 
of the norm transfer through the DB), it is sufficient to 
consider only two nonlinear sites with a third linear site 
attached. This is a strong indication that the destabi- 
lization threshold during the collision of the two objects 
can actually be linked to the Peierls-Nabarro barrier of a 
single DB. 



7. Possible Applications 

To start with, we would like to comment on the va- 
lidity of the DNLS Q to describe actual experiments of 
BECs in OLs. Experimental realizations have been per- 
formed for values of A = 2/7/ J in the range 10 5 -10 3 , 
while the number of atoms is typically N ~ (9(10 4 — 10 5 ). 
These estimations lead to experimentally feasible param- 
eters A = X/M < 1 for which the DNLS is a good ap- 
proximation. For example, the experiment of |8J shows 
that the the BEC dynamics in an OL with parame- 
ters N « 2 x 10 5 , J = 0.UE R , 2UN w 12E R (where 
Er = h 2 k\/(2m) is the recoil energy and Rl is the laser 
mode which traps the atoms), and M = 200 wells is de- 
scribed very well by the DNLS with effective nonlinearity 
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Figure 5: The destabilization of the DB in the trimer (solid line) is compared to the trimer with a third linear site (dashed 
line). On the timescale, where the destabilization takes place, the dynamics is qualitatively very similar (see enlargements in 
upper pictures). Total energy is E t — —1.31 > ^thrs and A = 3. (a) The initial condition (18) is given by 5 A — 0.016, 5$ = 7r, 
which is the same as in Fig. pp. (b) The initial condition is determined by the parameters 5 A — 0.152, 5$ = tt/2. 



A w 0.5. To estimate the minimum duration of an ex- 
periment with BECs probing the destabilization process, 
we rewrite our dimensionless time t — Jr in terms of 
real time r. Typical values for J are taken from |42| : 
The Josephson energy Ej/ks = 378 nK leads to the tun- 
neling rate J/H = Ej(Nh)- 1 = 16.5 Hz (for N = 3000 
atoms). Hence, with these parameters, t = 10 (which is 
a typical timescale after which the destabilization process 
took place) relates to r w 0.6 sec. 

The thresholded transfer of norm through a bright 
breather that we analyzed may lead to interesting appli- 
cations for blocking and filtering atom beams. It could be 
a powerful tool for controlling the transmission of matter 
waves in interferometry and quantum information pro- 
cesses [43J. 

In a similar way, our findings can be related to the 
field of optics, as the DNLS is capable of describing wave 
motion in nonlinear optical waveguide arrays. Discrete 
breathers in such two-dimensional networks were inves- 
tigated in the past years both theoretically and exper- 
imentally H3 OH l3^H35l [441 145] and can exhibit a rich 
variety of functional operations such as blocking, routing 
or logic functions [33j|34]. Experimental evidence of the 
destabilization process of the stationary DB should be 
observable in nonlinear waveguide arrays and might lead 
to functional operations such as filtering optical beams. 



Moreover, in view of a molecular trimer, applications 
in terms of targeted energy transfer (introduced in [JHl 
47J) with a threshold are conceivable, e.g. in the field of 
biophysics or biomolecular engineering. 

8. Conclusions 

The threshold and the tunneling process during the 
collision of a DB with a lattice excitation (e.g. a moving 
breather) was described analytically by defining the 2- 
dimensional Peierls-Nabarro energy landscape. The PN 
landscape restricts the dynamics of the trimer and the 
accessible region of the phase space. This restriction of 
the dynamics becomes very pronounced at the destabi- 
lization threshold, which is identified with a rim in the 
PN landscape. The effect is found for a broad range of 
the phase difference £</>i2 between the colliding objects. 
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